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Q: 

^ , A formula for an improved spin-dependent potential between a heavy quark 

^ ' and a heavy antiquark was developed using Heavy Quark Effective Theory tech- 

^\ , niques. The leading logarithmic quark mass terms emerging from loop contri- 

' butions were explicitly extracted and summed up. There is no renormalization 
scale dependence in this new formula and it includes both the Eichten-Feinberg 

■ formula as well as the one-loop QCD result as special cases. The hyperfine 

ly-^ I splittings were calculated using the spin-spin part of this improved formula. For 

^\ • charmonium the J/ip-rjc S'-wave splitting agrees well with the data while the 

I energy difference AMp between the center of gravity of the l'^Po.1,2 states and 

^ • the l^Pi state has the correct sign but is somewhat larger than the experimental 

I I value; however, there are also several other contributions to AMp of comparable 

magnitude (~ 1 MeV), which we discuss. The corresponding mass differences 

(— I I for the bb and be mesons are also predicted. 



Recently, significant progress has been made in the theoretical study of the spin- 
^ dependent potential between a heavy quark and a heavy antiquark. An improved 

5^ ! formula for the spin-dependent potential was developed 111 and used to calculate the 

hyperfine splittings in the cc, bb, and be systems i. The improved spin dependent 
potentials were derived from QCD first principles using the techniques of the Heavy 
Quark Effective Theory (HQET) i. The spin-dependent potential was separated into 
short distance parts involving Wilson coefficients and long distance parts which were 
expressed in terms of gauge invariant correlation functions of the color-electric and 
color-magnetic fields weighted by the Wilson loop path integral If the tree level 
values for the Wilson coefficients are used the potential reduces to Eichten's and Fein- 
berg's result i. And using the one-loop values of the Wilson coefficients, also calculat- 
ing the correlation functions to one-loop in perturbation theory, the spin-dependent 
potential at the one-loop level in perturbative QCD Si is recovered. However, the 
leading logarithmic terms appearing in perturbative calculations were also summed 

*Presented at the International Symposium on Heavy Flavor and Electroweak Theory, Beijing, 16-19 
August 1995. 



up in Ref. [1] using the Renormalization Group Equation (RGE) to obtain a scale 
independent result. Therefore, the spin-dependent heavy quark-antiquark potential 
derived in Ref. [1] is scale-indeDendent and thus improves upon and generalizes both 
Eichten's and Feinberg's result a and the one-loop perturbative result Q@. In addition, 
this improved result lil satisfies all the general relations among the different parts of 
the spin-dependent potential 0. 

The general formula for the spin-spin part of the renormalization-group-improved 
spin-dependent potential that was derived in Ref. [1] was used to calculate the hy- 
perfine spin splittings in the cc, bb, and cb systems in Ref. [2]. Since the spin-spin 
potential is a short distance feature, perturbation theory can reliably be used in the 
calculation. The result for the l^Si — I^Sq splitting between the J/tjj and the rjc 
agrees well with the experimental value i and the predictions for the mass differences 
ip' — rj'^, T{lS)—rib, T{2S)—ril, and B* — Bc are reasonable. However, the contribution 
to the P-wave energy difference, AMp, between the center of the gravity of l^Po,i,2 
states and the l^Pi state, while having the correct sign, is somewhat larger than 
the experimental data i . That is, when the contributions of the leading logarithmic 
terms are summed up and included, the agreement with that data is not as good as 
when only the one-loop perturbative spin-spin potential, in which the leading log- 
arithmic contributions are not summed up and included, is used. There are several 
other contributions to the rather small energy difference AMp (~ 1 MeV) which es- 
timates indicate are of the same order of magnitude as the spin-spin contribution. 
It therefore appears that the agreement of the one-loop perturbative result with the 
data is probably fortuitous. 

In the derivation the renormalized two-particle effective Lagrangian was first 
calculated to order Then, treating the terms of higher order in 1/m in the 

effective Lagrangian as perturbations, the four point Green's function on the Wilson 
loop i with the time interval T was calculated in the limit where m — > oo first 
followed by T — i> oo &0. In this limit, using standard perturbative methods, the 
large T behavior of the Green's function is of the form 

/ oc e-^^^'"'"). (1) 

From Eq. e{m,r), the potential energy between the quark and the antiquark, 
can be extracted. Expanding e(m, r) in powers of 1/m each of the spin-dependent 
potentials can be factorized into a short distance part, involving Wilson coefficients, 
and a long distance part, which can be expressed in terms of correlation functions 
of the color-electric and color-magnetic fields weighted by the Wilson-loop integral. 
Using the notation of Ref. [1], the resulting spin-spin potential is 

Si ■ S2 r 

AHss{mi,m2, r) = C3(/i, mi)c3(/x, m2)V4,{n, r) - 6Ncgs{fi)d{fi)6'^{r) (2) 

6mim2 L J 

where mi, m2, and Si, S2 are the masses and the spins of the heavy quark and the 
antiquark, respectively, fi is the renormalization subtraction point, is the number 



of colors, and 5's(/i) is the running coupling constant. The Wilson coefficients C3(/i, m) 
and d{fi) were calculated in leading logarithmic approximation in Ref. [11] and Ref. [1], 
respectively, and are 
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In Eq. (|^) V4(/i, r) is the color magnetic-magnetic correlation function which can be 
expressed as 



V4(/i,r) ; 
where (■ ■ ■) is defined by 
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Here C{r,T) represents the Wilson loop S, P denotes the path ordering, and r = 
|xi — X2I. 

We emphasize that this is a general result for the hyperfine part of the spin- 
dependent potential to order 1/m^. It absorbs the short distance contributions to the 
potential into the coefficients c^{^,m) and d{ii) while the long distance contributions 
to the potential are contained in the correlation function V4(/i,r). Moreover, the re- 
sult is independent of the factorization scale since the /i-dependence in the coefficients 
cancels the /^-dependence in the correlation function. The ffist term in the bracket 
in Eq. @) is a nonlocal term while the second term is a local one which is generated 
by mixing with the ffist (nonlocal) term under renormalization. We note that if the 
coefficients are evaluated at tree level- i.e., C3(/i, m) = 1 and = 0, the potential re- 
duces to the Eichten-Feinberg result And if these coefficients are expanded to order 
as{^) and the correlation function is also evaluated only to one-loop, the logarithmic 
terms in Eq. (H) then reduce to the one-loop spin-spin potential M. Therefore, this 
renormalizat ion-group improved potential, Eq. (2), extends both Eichten's and Fein- 
berg's result as well as the one-loop perturbative potential, containing each of these 
results as special cases. 

To ffist order perturbation theory in AHss the energy shift caused by AHss{r) is 
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dh^l,(r)AHJr)^,,(r) 
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where \l//^/^(r) is the nonrelativistic wavefunction of the bound state with total angular 
momentum I and 2;-component For simplicity we suppress spin and color indices 
and retain only the space-dependent indices. Separating the radial part, u{r), we 
write \E'«,i,(r) as 

vl/,,,^(r) = w(r)y,,je,0), (8) 

where Yi^rn{0, 0) are the standard spherical harmonics. 

The radial wavefunction was obtained by numerically solving the Schrodinger 
equation. For comparison, we used three popular potential models. One was the 
Cornell model 111, the second one was the logarithmic potential 111, and the third one 
was the improved QCD-motivated potential 0. 

The calculation also required an expression for the running coupling constant 
as{q). The familiar RGE, one-loop result is 

(^s{q) = (9) 

6oln-r2- 

MS 

where bo = llNc — 2Nf and Nf is the number of quark flavors. It is clear from 
Eq. (P) that as{q) contains a Landau singularity in the nonperturbative region when 
= A|j^ and becomes negative for < Afj^-. To avoid the resulting numerical 

iVL o iVi o 

ambiguities we first moved this singularity to = and used a modified form of 
as{q) in the actual numerical calculations; namely, 

Att 

as{q) = r-- (^'^) 

6oln -^ + 1 

\ MS / 

The value of Ajj^ was taken to be 200 MeV and 250 MeV in the numerical calculations, 
which is within the experimental range = 195 + 65 — 50 MeV i. 

To explore the sensitivity of our results to the location of the Landau singularity 
in as{q) we replaced the expression, Eq. (|TUp, by 

Att 

as{q) = ^ (11) 

V MS / 

and varied A^. The results for the 5*- wave hyperfine splitting were not sensitive to 
and only for large A^ did AMp significantly decrease. To fit AMp to the measured 

value required A^ quite large, about 16, clearly out of the perturbative region. 

Our numerical results for the ^5*1 —'^ So and ^Pj—^ Pi splittings, where ^Pj denotes 

the center of gravity of the '^-Po,i,2 states, for the three potentialgiHlEJ'EJ are presented 

in Tables I, II, III, respectively. For comparison we have also included the results for 

the 25* and 2P states. 



Table 1. The hyperfine spin splittings in MeV predicted using the Cornell potential 



cc bb be 

Aps (MeV) 200 250 200 250 200 250 

E{l-^Si) - E{1^Sq) TTTT 128!3 97?7 1043 67l) TLI 

E{2^Si) - E{2^So) 75.3 82.5 39.6 42.2 37.7 40.3 

E{l^Pj) - E{l^Pi) -4.8 -7.0 -3.0 -4.1 -4.0 -5.6 

E{2^Pj) - E{2^Pi) -3.6 -5.2 -2.1 -2.9 -3.0 -4.2 



Table 2. The hyperfine spin splittings in MeV predicted using the Logarithmic potential 

|13 

CC bb be 

Aps (MeV) 200 250 200 250 200 250 

Eil'^Si) -E{l^So) 106.1 117.1 36l) 38^2 4(16 42!6 

E{2^Si) - E{2^So) 54.2 59.7 18.7 19.8 21.2 22.3 

E{l^Pj) - E{l^Pi) -5.4 -7.8 -3.4 -4.5 -4.5 -6.4 

E{2^Pj) - E{2^Pi) -2.6 -3.8 -2.1 -2.8 -2.8 -4.0 



Table,3. The hyperfine spin splittings in MeV predicted using the Improved-QCD motivated poten- 
tialllil 





ce 




bb 




be 




Ams (MeV) 


200 


250 


200 


250 


200 


250 


Eil^Si)-Eil'So) 


107.9 


119.1 


44.6 


47.6 


43.4 


45.7 


E{2^Si)-E{2^So) 


68.5 


75.6 


20.9 


22.4 


25.2 


26.7 


E(l^Pj)-E{l'Pi) 


-4.6 


-6.7 


-2.7 


-3.7 


-3.7 


-5.3 


E{2^Pj)-E{2'Pi) 


-3.4 


-5.0 


-1.8 


-2.5 


-2.6 


-3.8 



The main features of these results can be summarized as follows: 

• The results are /^-independent, as they must be. 

• The calculated energy difference between the J/"^ and the rjc mesons is quite 
close to the experimental value i for all three potentials. 

• For each of these three potentials we predict the energy difference between 
and 7]'^ to lie within the range 55 — 80 MeV. 

• For the bb system there are significant discrepancies between the Cornell model 
and the other two models for the S-states. Since the Cornell model, with the 
parameters fit to the cc spectrum, does not predict the bb spectrum very well, 
the results calculated in the other two models are probably better predictions 
for the energy difference between the T(IS') and the rib (35 — 50 MeV) and 
between the T{2S) and the r/^ (20 MeV). 



• The predicted energy difference between B* and Be meson is in the range 40 — 70 
MeV from all three of these models, which is consistent with previous results 0. 

• The calculated value of AMp = E{l^Pj) — E{l^Pi) for the charmonium IP 
states is in the range of —4 to —6 MeV, which has the same sign but is sev- 
eral times larger than the experimental value of —0.9 ± 0.2 MeV i. This is 
not surprising since there are several other contributions to AMp which esti- 
mates indicate are comparable in magnitude to the contribution coming from 
the hyperfine spin-spin interaction, Hgs- In fact, it is surprising that the pre- 
diction from only the one-loop spin-dependent potential is quite close to the 
experimental data. 

To summarize, the hyperfine spin splittings in the cc, bb, and be system were cal- 
culated i using the RGE improved perturbative spin-spin potential i. The results for 
the hyperfine splittings of the 5- wave states agree with the J/^ — r/c measured split- 
ting u and the prediction for T — r]b splitting is reasonable. However, the contribution 
to AMp = Ei^Pj) — E{^Pi)foT the charmonium P-wave states is somewhat larger 
than the experimental data eI, although it agrees in sign. That is, after summing up 
the leading logarithmic terms and including them in the perturbation calculations, 
the agreement with the data is not as good as the one-loop calculations QBOIli. But 
there are several additional contributions that are possibly comparable in magnitude. 
These include the following: 

• The contributions of the spin-orbit and and tensor potentials in the second 
order of perturbation theory: These contributions to AMp only cancel to first 
order in perturbation theory. However, according to the power counting rules 0, 
the spin-orbit and tensor potential potentials shift the energies of the P-wave 
states by an amount of order mv'^ in first order, which indeed cancel in AMp, 
but they do make a contribution to AMp of order mv^ in the second order of 
perturbation theory. This estimate is several MeV for the P-wave charmonium 
states, and therefore should not be ignored. 

• Higher dimensional operators: Unlike the dimension-six operators, these give 
non-zero contributions to AMp even at tree level. Compared to the one-loop 
contribution, these are suppressed by but enhanced by and v^/as ~ 1 in 
charmonium. 

• The color-octet S-wave component in P-wave quarkonia states 111: This com- 
ponent of the wavefunction receives a tree-level contribution from the local 
term 6^{r) in the spin-spin potential. This contribution too could be of order 
v'^/ag ~ 1 compared to what has been calculated. 



• Next-to-leading order contributions from the two-loop potential: These are 
suppressed by order a^, but since is not a very small quantity in charmonium, 
one cannot dismiss the possibility that this contribution could be significant. 

Before comparing with the experimental value of AMp in charmonium, which is only 
about 1 MeV, all the above contributions should be included since they are possibly 
comparable in magnitude. In the bb case these effects are less important and one can 
expect the perturbative calculations the bb system to be more reliable, although AMp 
is smaller, also. 
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